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I. INTRODUCTION 


- nancies description of the interaction of oe weet matter 
over a large range of energies is a eipeees of basic importance i in a wide variety | 
of theoretical and applied areas. From the theoretical standpoint, aids Gen 
of energy loss and range of electrons in many different materials have formed 
the basis of at least two extensive tabulations. me Both of these works are 
restricted to electron energies 2 10 keV and are based on the Bethe theory of 
stopping power including various modifications and corrections (e.g. density- 
effect corrections). We feel that similar tabulations for electron energies 
< 10 keV, based ona priori calculations using currently available thearetical 
information, will provide useful guides far interpretation of experimental data 
as well as input for calculations in applied areas. 

Our work here will involve model calculations to describe the extended 
electron states of a solid (valence band or conduction band). The more tightly 
bound, inner shells of the atoms in the solid will be assumed to be essentially 
unchanged in character from those in free atoms. Thus, excitation of electrons 
from the inner shells will be based on calculations of atomic, generalized 
oscillator strengths (GOS's). More eeveasny: the differential inverse mean 
free path (DIMFP), which forms the basic function required in our work, will be 
derived from an electron gas model descrihing the conduction band electrons in Al 
and from a model insulator theory applied to the valence band in Al,0,- The 


DIMFPFP's for electron interaction with the inner shells of the Al and O atoms 


lehenti ection bite khan ch-3 
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are calculated from a priori atomic GOS calculations for exciting electrons to 
the contiruum from the 2s and 2p sub-shells of aL and the K shell of o.4 The | 

DIMFP far removing an electron from the K shell of Al is obtained from GOS 

values based on hydrogenic wanetuncGons.” Given the DIMFP's associated with 

the most important electron interaction processes in the solids we then calcu- 

late inverse mean free paths, stopping powers, csda ranges, and range ard en- 

ergy straggling fcr electron energies from a few eV to 10 keV. 

The details of the components of our calculations are described more 
fully in the next four sections. In section VI exchange corrections are dis- 
cussed, expressions given for the exchange corrected DIMFP's, and formulas 


used in these tabulations are displayed. 


i. GENERAI, FORMULATIONS 


A charged particle passing through a solid interacts with a large aumber 
of electrons simultaneously and it is thus appropriate to speak of a mean free i 
path of the charged pabtdete against energy loss to the solid. Assuming the ef- 
fect of the charged particle on the medium may be treated in first Born approx- 
mation, the inverse mean free path, differential in momentum transfer, hk 5 


and energy transfer, hw, for a particle of velocity v is given by 


ai ze zy ] 
se = Im} ——- (1) 
daw ake €(k, W) 


where €(k,W) is the exact dielectric function of the saat We assume in 


this work that the solid is isotropic and homogeneous. 
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For our calculations of inverse mean free path, s topping power, etc., © 
it is sufficient to compute inverse mean free paths differential in energy 
transfer only. This differential inverse mean free path (DIMFP) for energy 


loss hw by an electron with energy Es mere in the solid is given by _— 


. k 
om 2 tf * & | Jia 
T(E, hw) ana Ta E k Im ck (2) 


where hk | =V2m [VE + VE-he| anda, = ohne This expression 
assumes that the energy-momentum relation far a swift electron in the solid 
does not differ appreciably from that of a free electron in vacuum. 

Given ek,w) for the solid, the quantities of interest here follow directly 
from T(E,hw). The inverse mean free path of the electron, p, is given by 


integrating over allowed energy transfers as 


BK(E) = f dhw) T(E, nw). (3} 
The rate of energy loss of the electron, ar the stopping power of the medium, 
is given by 

SE) = -dE/dx = fd(hwi hw T(E,hw) , 


and the mean square energy loss per unit path length by 


27 (BE) wf acneyenes? T(E, hw). (5) 


With these results we may calculate the range of an -lectron in the continuous- 


slowing-down approximation (csda range) by 


ee ee Bed 


wi ee 
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E 
RAE) =f aE’ /SE') . | (6) 


E 
3 


The lower limit on this integration will be discussed further in Section VI. The 
mean square fluctuation in the range or "range straggling" will be calculated 


from Eq. (3) and Eq. (4) as® 


E 
° 


2 2 3 
(R-R) =f dE AN(E/ CSE)” (7) 
AY 
In practice, the DIMFP will be evaluated as a sum of contributions from 
various distinct processes. For example we calculate a DIMFP for removing an 
electron from a given inner shell, a DIMFP far plasmon excitation, etc. The 


total DIMFP used to describe the interaction of an electron with the given solid 


will ke given by 


TE,W) = 2 7,(E,«) (8) 


i i 
where the sum over i adds the contributions from the various interaction proc- 
esses. The evaluation of the 7's far Al and Al,0, is described in the next three 


sections. 


Il. DIMFP's FOR THE ALUMINUM CONDUCTION BAND 


The conduction band of Al will be described by an electron gas model. The 


dielectric response functicn in the form given by pinihara” is 


e(k, a) = 14 (x°/2") [LE (x,2) + if, (2) ] (9) 


¢ 


in terms of the dimensionless variables x = hw/E, and z =k/2, withE, the 
Fermi energy ard k, the Fermi wavenumber of the electron gas. Also, 

2 2 : . ‘ So ae 
xX = e /whv, where v, is the Fermi velocity. The function f is given by 


_ . a, : 2. ors ae _ 
£ (xz) = & + = {1 -(z - x/4z) HN eens a | 


A 7 2 z+x/4z41 
+ at} (z + x/4z)" Jén | Saar Ls | (16a) 


The function f, is defined by 


x/8z , for z+x/4z < 1 
£,(%2) = ZL - (2 -/42)71, for |z-x/4z|<1<z+x/4z (10b) 
0, for |z-x/4z/>1 . 


From Eg. (2), the DIMFP is given by 


2 + zf, 
T(E Buh = Fe =f Sei ge ie ee a as. 
Ore. (z +x £.) +(x f,) 


where z, = ky /2k, . Energy losses of an incident electron to an electron gas 


may be divided inte: (a) losses resulting from excitation of single electrons out 
of the Fermi sea and (b) losses to collective oscillations of the electron gas 
(plasmon excitation). The contribution to the total DIMFP due to (a) is identi- 


fied with the region in which f, is non-zero. This electron-electron term is 


written as 


S(/xtl +1) 
duos 4 ; : ? a e- x a zf, sate 
dx Le Gas aid Ta (e+ 1) z 


2 2 
dori -1) (tx A) HOE) 


4 
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where € is the incident electron energy measured from the Fermi level in 


Ca erre 


units of the Fer:ni energy, i.e., ¢ =(E-E,)/Er, and the step function — z 
@(€ - x) restricts the incident electron to states above the Ferrai level. ! 

The plasmon contr bution arises from integration of Eq. (11) in the 
region where f, > 0. The integral is zero except on the plasma resonance 


line defined by z = z6) through the equation 
- 2 2 £ : 
¥(x,2,) = 2. +x £,(%25) = 0. (13) : 


The DIMFP for plasmon excitation is found to be : 


z 
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where x. and x are the solutions of the equations 
min max 


eee ele 


Xin 7 42 6(%min) [wetl - 2 5(%min) (15) 
and 
Mer I {1+ Zax) , (16) 


ee 


For our later discussion of exchange carrections we need the DIMFP far 


creation of secondary electrons in this model. The results of Ref. (7) lead to 


min(¢’ +1, €) Zz, 


2 - 

x dz 3 
T(E,€) = fx f (17) 7 
= Scere e’ z- 2* | €(2k,2,E,x) iy . 


where minfa,b) represents the smaller of the quantities a and b. 
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The use of this electron gas model, though based implicitly on the | 
assumption of high electron densities, gives unexpectedly good results for 


real metal conduction band densities. In addition, Eq..(1) is derived in first 
i 


Born approximation and would be expected to fail when the velocity of ia ies 
cident alectron approaches that of electrons in the electron gas. Lindhara, 6 : 
however, points out that this equation may be reasonably good at axy velocity 
since: (a) the relative velocity of the incident eiectron and gp clear ba 
wiacivonst in the elec izon gas remains of the order of the Fermi velocity no! 
matter how slowly the incident electron moves, and (b) at iow velocities “ 
Coulomb field of the incident electron is strongly screened by polarization in 
the electron gas so *1t for many purposes it can be considered smal? in the 


sense of perturbation theory. 


IV. DIMFP'S FOR INNER SHELLS 


From a general expression for the dielectric function of a homogeneous, 


t 


i 5 
isotropic egetan we may show for values of w which correspond to ionization 


of the a inner shell in a solid that 
5 
-] 2mn.e df (k,w) 

im | Sow] ~ Im €(k, W) = Tae. apg (18) 
where dé /du is the generalized oscillator strength (GOS) far transitions from 
the gh independent inner shell. Here n, is the number of a inner shells per 


unit volume in the given solid. Equation (2) thus leads to 


wa Barta i aca aor reacia ve 


ll 


T(E, hw) 


81a | *n, df (k, ) 
al 3 Ci 2 


ERK have a( ha) 


where hw is the energy transfer and R =e ?/2a 

Generalized oscillator strengths for the ionization of L sheil electrons in 
Al have been calculated by Manson’ using a nonrelaSivistic Hartree-Slater central 
field model of the atom. These GOS values have been used as input for numerical 
evaluation of the integral over momentum transfer in Eq. (19) to obtain differen- 
tial cross sections, T/ n= do. /A( hw), for the 2s and 2p subshells of Al. These 


cross sections when multiplied by the appropriate value of n. give the DIMFP for 
Al, Al,O 


2s would denote the DIMFP for 


the particular material. For example, T 
removing a 2s electron from an Al atom with the density of Al atoms correspond- 
ing to that in Al,0,.- This same type of nctation will be used for contributions tc 
the inverse mean free path and stopping power. The binding energy for 2s and 2p 
electrons in Al is takentobe 118.524eVand 80.88 eV, respectively, for 
these calculations. 

Similar calculations have been done for the ls shell in oxygen. The values 
of the GOS for excitation to the continuum from the Is shell of Owere taken 
from the work of McGuire. * The binding energy of the 1s shell in O is taken to 
be 536.6 eV. 


Our calculation of the DIMFP for excitation of electrons from the K shell 


of Al to the continuum is based on cross sections derived using hydrogenic wave 


functions. A detailed discussion of this type of cross section calcuation is 
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given in Ref. 5. We have used the equations in this reference to calculate the 
DLIMFP based on a binding energy far the K shell in Al of 1545.8 eV and an 


"effective atomic number" of Z. = 12.5643. 


. V. DIMFP'S FOR VALENCE ELECTRONS IN A1,0, 


Leyes 


Since A1,0, is a good insulator with a band gap of about 9 ev'0 it is 
desirable to represent the response of its valence electrons on the basis of 
a quite different model than that used for the conduction band in Al. The 
model which we have developed for this purpose is related to that employed 
by Fry! in which the ground state wave function of the valence electrons is 
described in the tight-binding approximation, while excited states are repre- 
sented by orthogonalized plane waves (OPW). In our use of the model to ob- 
tain a dielectric response function we fix the normalization of the OPW excited 
states by requiring that the sum rule f dw wimf[ €(k,w)] = 2m ne-/m . 
where n is the density of electrons in he valence band. In addition we assume 
that the solid is unifarm ;: nd homogeneous. The dielectric response function 
corresponding to this model solid is corvenient and flexible for use, can be 
fitted to the optical dielectric function in the limit of very long wavelengths 
(k°0), and describes the single-particle properties of excited electrons. The 
existence of plasma oscillations emerges naturally as one studies the response 
of the system to longitudinal electric perturbations. A 


Since a detailed discussion of the insulator model is planned for publi- 


eon we quote here only the results needed for these calculations. The 
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result required here is the imaginary part of the dielectric response function 


for the model insulator given by 


Im[e(k,w)] = mmeT /hBkA (20) 
| where 
rei a | 
{a+ (k-p) } {a +(kt+p) } 
_ 32a° [ 1 fe 1 | (21) 
(is py (aa ey eo 4(k=p)" a etkani 
, —-lo24pk a 
(vey Gxeicy: 
np,2, 27f 1 32a" 
A = [w+ a (k “+a i) Soo a a z| (22) 
(40° +k ) 
and 
1/2 
p = [m(w-ayVh8 J : (23) 
Here hw, is the excitation energy of the valerce electrons and B ard a are 


| parameters which may be adjusted to make the theory agree with optical dielectric 


function measurements in the k°O limit. In the k>0 limit we have 


[ 2" wne- on p° 

Im € oO, W ee ree ; 24 

( ] 3np i. 2B a’) (a 4p) (24) 
5 m 


Given the imaginary part of the dielectric function, Eq. (20), for fixed values of 


n, B, Ws, and @ the real part of €&k,W) may be obtained mmerically using the 
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Kramers-Kronig relation, 


sad W Im[€(k,w)] 
Ref €(k,w)] = 1+ 2 f dw, _—— = : (25) 
o (ey 


° 


_ We fit Eq. (24) to experimental measurements of the imagivany part of 
the optical dielectric function as measured by Arakawa and Williams m with 15 
valence electrons per Al,O, unit, density of Al,O, of 4.05 g/m’, hws, = 9eV, 
B = %, and aa = 0.78. Given these values, Im[ - a ] is calculated 
from the results of Eq. (20) and Eq. (25). The remaining 9 of the 24 valence 
electrons in Al,O, are taken to form a second tight-binding level with hug = 


29 eV, B=%, and aa, = 1.6. 
VI. EXCHANGE CORRECTED DIMFP'S AND FORMULAE FOR THE TABULATIONS 


We have included the effect of electron exchange in our calculations ina 
simple manner based on the form of the Mott formula (non-relativistic Mdller 
formula) for scattering of an incident electron with a free electron. The cross 


section for finding a scattered electron with energy in the interval W:W+dW is 


given by” 


4 
#2 . se / 41-1 | si 
dw E Ww (E- W) W(E-W) 


for an incident electron of energy E, except for energies close to W = 0 and 


W=E. Near W =0 and W = E the interference term (third term on the right side 


of Eq. (26))is effectively 7-_: o. 


15 
The DIMFP for excitation of an electron from a particular state i may 


be weitten in the form 


1 
7 Fi(E, hw) . (27) 


T(E, hw) = 
If we assume that the width of the level from which an electron is excited is 
quite narrow, we ov‘-in from Eq. (27) the DIMFP for production of a second- 


ary electron with energy E, as 
P(E,E.) = + F.(E, EX¢+E) (28) 
gee OE ge Ss 


where = is the binding energy of the = state (a positive quantity). The ex- 


change corrected DIMFP is taken as 


exc m5 eaee 
7; (E, hw) = E |r ce, nw) + FE, E+E: hw) 


1/2 
- [2 NEFF] eeneye,ce, E+ E; - na] + (29) 


Since ET, = 1/(hw)* for large E and hw, Eq. (29) reduces in this limit to the 
form given by Eq. (26). The factor 1 -NES 7E reduces the contribution of the 
third term in Eq. (29) as E> ES - This form far the exchange corrected DIMFP 
has been used in our calculations for all the inner shells and for the two valence 
levels in Al,0, (since our model assumes the width of these levels to be quite 
na~row). 

If we now define the more energetic of the two electrons after collision 


to be the primary, and account for exchange through Eq. (29), Eq. (3) gives the 


LAR ea oe ake ie aS 
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the contribution to the inverse mean free path due to excitation of an electron 


from the ith level as 


(E+E*)/2 
HAE) .| d(hw) 7. (E, hw) : ve 20% (30): 
Et fe tae : 
Similarly for the stopping power and mean square energy loss per unit path 


length we have from Eq. (4) and Eq. (5) 


(E+E.)/2 

S(E) = a(hw) hw 7, (E, bw) (31) 
8 ‘ 

Ei 
and 

(E+ E?)/2 

07 (E) = a nay(nw)” 7. (E, nw) ' (32) 
ES 
2 


The exchange correction for excitation of electrons from the Al conduction 


band of finite width requires a slightly different form. We take 


-1/2 v2 


tl e,x) = Tae (ErX) + T.(€,€-x) - [1 - (et) IL tr fexT (6 e-x) J (33) 


ee 


where ee and fet are defined in Eq. (12) and Eq. (17). Thus the contributions to 
the inverse mean free, stopping power, and mean square energy loss per unit path 
length due to excitation of electrons from the conduction bard in Al are calculated 


from 


exc 
Bl ©) -f dx TI, (6x) » (34) 


AAS RS A Serene RRS BP ae AER hl TN i AF OA RAEN leet eee senemnnmanmmemamnenie en a owen 
: 3 eT re re” a 
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x/2 
s_(e) = ex x ei (ex), (35) 
fe} 
ard 
x/2 
a tee dx x2 - se (€,x) . (36) 
ee ee 
fe} 


No exchange correction is applied to To Eq. (14). The inverse mean free 


path contribution due to plasmon excitation is given by 


x 
max 


pyle) =f axe Tylerx - (37) 


x. 
rain 


The contribution to stopping power and mean square energy loss per unit path 


length are given by 


x 
max 
$i°) .{ clear aa a (38) 
x. 
‘min 
and 
x 
max 
a2 (ec = d ar (€,x) 39 
pl ) ve x pl > ad ( ) 
*min - 


For the remaining calculations we form the sums 


Sexc(E) = % 8,(E) (40) 


Ye LN a SE PES 5S 


~ Spl de OF Ae 


a | 18 


and 


a a® (£) - Da? ey . (41) 
1 


where the index ‘ includes the terms appropriate for a given solid, including | 
exchange corrections as indicated above. The csda range is calculated from, 
: R19) F) -f dE‘/S_(E) | (42) 
10 eV 


carresponding to an electron slowing down in a continuous manner from an energy 
E to 10 eV. The mean square fluctuation in the csda range based on Eq. (7) is 


calculated as 


E 
[AR 10) a .| dE’ @ EV IS (EVP , (43) 


10 eV 
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Vl. ALUMINUM: EXPLANATION OF TABLES 


GENERAL NOTES 
1. Electron energies are measured from the top of the conduction 
bard (Fermi ee The Fermi energy is taken tobe E, = 11.6 ev. 
2. The density of solid Al is taken to be 2.71 heta : , 


3. The computer-printed units are translated as: 


¢ EBV sed e0sierwes eV Ass evreence cece A 
EV2 hawucces (ev)? 1 es ene, a: 

3 2 

G/CM3..... ecm Az: sess A 


4. The numerical print-c:1t is in the form, e.g. 


2.8D-01 = 2.8x10. 


TABLE 1A - INVERSE MEAN FREE PATH OF ELECTRONS IN ALUMINUM 


EL ~ EL Hoe as given by Eq. (34) 
PLASMON Bot as given by Eq. (37) 
Al ‘ ‘ 
AL(n4£) Hig as given by Eq. (30) with Eq. (19) 
INVERSE MFP p. ~ total inverse mean free 
path = sum of individual 
contributions. 


TABLE 1B - STOPPING POWER OF ALUMINUM FOR ELECTRONS 


EL-EL Sue as given by Eq. (35) 


21 


‘PLASMON . S.4 as given by Eq. (38) 
fl : ave 
AL(n2Z) nt as given by Eq. (31) with Eq. (19) 


STOPPING POWER S$ - total stopping power = sum of 
; individual contributions 


TABLE 1C - CSDA RANGES AND STRAGGLING OF ELECTRONS IN ALUMINUM 


CSDA RANGE Rio) 7 the range of an electron in the 

(E TOW EV) cbttdmious -slowing-down approximation 
in going from an energy E to 10 eV, as 
given by Eg. (42). 


STRAGGLING 


MEAN SQUARE a ~ the mean square fluctuation in 
ENERGY LOSS the energy loss per unit path length, as 
given by Eqs. (41), (36), (39) and (5) 
plus (19). 
MEAN SQUARE [AR ign - the mean square fluctu~ i 
RANGE FLUCTUATION atioh in the range about the mean csda | 
range R, 10)’ as given by Eq. (43). 
| 
2431/2 
RELATIVE RANGE ALAR yo) Jat Rao) | 
} 
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IX. ALUMINUM OXIDE: EXPLANATION OF TABLES 


GENERAL NOTES: 


1. Electron energies are measured from the bottom of the conduction 


band. 


2. The density of solid ALO, is taken to be 4.05 gjen. 


3. The computer-printed units are translated as: 


AL203 ..... AlO, 


-2 
4. The mmerical print-out is in the form, e.g., 2.8D-2 = 2.8xl0 ~. 


TABLE 2A - INVERSE MEAN FREE PATH OF ELECTRONS IN Al,O, 


VAL (9 EV) 
VAL (29 EV) 
AL(nt) 


o(1s) 


INVERSE MFP 


Hg eV) as given by Eq. (30) +(2) +(20) +(25) 


M29 ev) as given by Eq. (30) +(2) +(20) +(25) 


ees as given by Eq. (30) +(19) 


Oiaide a as given by Eg. (30) +(19) 


p - total inverse mean free path = sum 
of individual contributions. 
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TABLE 2B - STOPPING POWER OF ALUMINUM OXIDE FOR ELECTRONS — 


TABLE 2C - CSDA RANGE AND STRAGGLING OF ELECTRONS IN Al,O 


VAL (9 EV) Sa ayy 28 given by Ea. (31) +(2) +120) +(25) 
VAL (29 ei 2 Sag ey 28 given by Eq. (31) +(2) +(20) +(25) 
AL(nt) ears as given by = (31) +(19) 
O(1s) go 293 as given by Eq. (31) +19) 


STOPPING POWER S - total stopping power = sum of 
individual contributions. 


3 


CSDA RANGE R - the range of an electron in the 

(E TO D EV) ch ous-slowing-down approximation 
in going from an energy E tc 10 eV, as 
given by Eq. (42). 


MEAN SQUARE a - = theimean square fluctuation in 


ENERGY LOSS the energy loss per unit path length, as 
given by Eqs. /41), (32) +(19), and (32) + 
(7) +{20) +(25). 


2 
MEAN SQUARE C4R 10 Ji, - the mean square fluctuation 
RANGE FLUCTUATION in the range about the mean csda range 
Roy’ as given by Eq. (43). 
2,1/2 
L A j 
RELATIVE RANGE {C Roy 4 i /R 9) 


STRAGGLING 


PATH OF ELECTRONS IN AL203 (DENSITY 4.05G/CM3) 
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